Fluctuating annihilation cross sections and the generation of density perturbations 
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Fluctuations in the mass and decay rate of a heavy particle which for some period dominates 
the energy density of the universe are known to lead to adiabatic density perturbations. We show 
that generically the annihilation cross section of the same particle also receives fluctuations, which 
leads to entropy perturbations at freeze-out. If the particle comes to dominate the energy density 
of the universe and subsequently decays, this leads to an additional source of adiabatic density 
perturbations. On the other hand, non-adiabatic density perturbations result when the particle 
does not decay but contributes to the observed dark matter. 
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I. INTRODUCTION 



Measurements of the cosmic microwave background ra- 
diation have revealed a highly uniform energy den- 
sity background with super- horizon perturbations on the 
order of one part in 10 5 . In the standard inflationary 
paradigm 0, > these density perturbations were created 
in the inflationary epoch when quantum fluctuations of 
the inflaton field expanded beyond the Hubble radius and 
were converted into density perturbations upon inflaton 
decay. However, to obtain the observed level of density 
perturbations from this mechanism requires tight con- 
straints on the inflaton potential 0. 

Recently, Dvali, Gruzinov, Zaldarriaga and indepen- 
dently Kofman (DGZK) proposed a new mechanism p| 
for producing density perturbations. A nice feature of 
their scenario is that the only requirements on the infla- 
ton potential are to produce the required e-foldings of in- 
flation and at a scale consistent with WMAP data. The 
DGZK mechanism posits the existence of some heavy 
particle S with a mass and decay rate that depend on 
the vacuum expectation value of some light field \. Here 
X is presumed to have acquired super-horizon fluctua- 
tions during the inflationary epoch; however \ never con- 
tributes significantly to the energy density of the uni- 
verse 1 . Nevertheless, the fluctuations in \ persist and 
result in fluctuations in the mass and decay rate of S, 
so long as the x mass m x is less than the Hubble rate 
H at the time at which fluctuations are transferred to 
radiation. In the DGZK mechanism the field S comes to 
dominate the energy density of the universe and decays 
into radiation while m x < H. Fluctuations in the mass 
and decay rate of S result in fluctuations in the duration 
of S energy domination, which in turn lead to adiabatic 
density perturbations since the energy of a massive S 
field redshifts more slowly than that of radiation. 

The DGZK mechanism has been studied extensively. 
For example, the evolution of the density perturbations 
that result from this mechanism has been studied in de- 



tail using gauge invariant formalisms in : TJ. These per- 
turbations are shown to possess a highly scale invariant 
spectrum in |8| and are shown to contain significant non- 
Gaussianities in The original DGZK mechanism has 
also been extended to apply to preheating as studied in 
|lfjj | . For discussions of the limitations of this mechanism 
see for example [TT| . 

In the original DGZK scenario 5] it is assumed that S 
decouples while being relativistic. In this paper we gen- 
eralize this to apply to the case where S freezes-out of 
equilibrium with a fluctuating annihilation rate (av) . We 
use the term "freeze-out" to refer specifically to the sce- 
nario where S decouples from thermal equilibrium after 
it has become non-relativistic. In this case the number 
density of S at a temperature T after freeze-out is 
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ns 
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(1) 



where m p \ is the Planck mass and ms is the mass of S. 
Therefore we expect fluctuations in the mass and annihi- 
lation rate of S during freeze-out to result in fluctuations 
in the number density of S. If S lives long enough to 
dominate the energy density of the universe and sub- 
sequently decays, these entropy perturbations are con- 
verted into adiabatic perturbations. These add to the 
ones produced by the original DGZK mechanism and the 
quantum fluctuations of the inflaton. 

This paper is organized as follows. In Section [H] we 
describe the density perturbations produced by our gen- 
eralized DGZK mechanism. Sections lllll and HVI contain 
explicit models for implementing our mechanism and for 
producing the fluctuating masses and coupling constants, 
respectively. Conclusions are given in Section In 
Appendix ^ an alternate analytical description is given 
which allows to track the evolution of the perturbations, 
while in Appendix [5] Boltzmann equations are derived 
and solved numerically to confirm the analytical argu- 
ments presented in other sections of this paper. 



II. ANALYTICAL DETERMINATION OF THE 
PERTURBATIONS 



1 The scenario where \ contributes significantly toward the energy 
density is called the "curvaton" scenario and was first proposed 



Our generalized DGZK mechanism includes a heavy 
particle S with mass ms, decay rate T and annihilation 
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cross section (cry) , where S decays to and interacts with 
radiation. We begin by identifying several key tempera- 
ture scales. The temperature at which S begins to ther- 
malize with radiation is denoted as Tthcrm- We assume 
for simplicity that S particles are produced only as they 
thermalize from radiation annihilation below T = Tthcrm- 
We also define: 

1. Tf.o.: Temperature at which S freezes-out of ther- 
mal equilibrium; 

2. 7d om : Temperature at which S begins to dominate 
the energy density of the universe; 

3. T(j cc : Temperature at which S decays. 

Since the number density of S particles falls off exponen- 
tially after S becomes non-relativistic, Tf. Q . is typically 
within an order of magnitude of 7715. Therefore in this 
paper we always take Tf. . ~ m<?. In terms of mj, T and 
(av) we also find 



Tic 
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where we have assumed Tdcc < Td om in the last equation. 
This condition is necessary for significant density pertur- 
bations to be produced by this mechanism. In Eq. (J2J 
the cross section is to be evaluated at the freeze-out tem- 
perature Tf. .. Note that for S particles to be produced 
in the first place we require Tthcrm > Tf.o.. 

As described in the period of S domination be- 
tween Tdom and Tdec gives rise to an enhancement of the 
resulting energy density compared to a scenario where 
the S domination is absent. Comparing energy densities 
at common scale factor one finds that after S decays 
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where 
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(4) 



and p ra d is the energy density which would result without 
any period of matter domination. As discussed in detail 
in Section IV, couplings to an additional field \ can give 
rise to fluctuations in ms, I\ and (av): 



m s 
(av) 

r 



m s (1 + 5 m ) 
(av) (l + 5 {av) ] 

T(l + S r ) , 



(5) 



where the barred quantities refer to background values. 
According to Eqs. 12I5J1 . these fluctuations give rise to 
fluctuations in Td om and Td oc which result in energy den- 
sity perturbations 



Sp 
P 



(6) 



Note that although Sp/ p contains no explicit dependence 
on S m , both 5i av \ and St are in general functions of 5 m . 

Comparing the energy density at a common scale fac- 
tor corresponds to choosing a gauge where the pertur- 
bation in the scale factor vanishes, ip = 0. Thus the 
fluctuation in the energy density computed here can be 
directly related to the gauge invariant Bardeen parame- 
ter [Tif 



C = -V- 

Thus we find after S decays 
1 „ 



Sp 
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(8) 



We can obtain the same result in synchronous gauge, 
where different regions all have the same global time. 
Since p ~ 1/t 2 in both matter and radiation dominated 
universes, one finds that Sp = on surfaces of constant 
time. Thus the Bardeen parameter is 



c 



Sa 
a 



(9) 



To obtain £, we only need to determine a(t, T, (av),m) 
and then compare two regions at fixed t, but different T, 
(av) and ms- Assuming the S particles freeze-out while 
non-relativistic and decay after dominating the energy 
density of the universe, this gives 

a(t) a(t dcc ) a(t dom ) a(tf. .) 



a(t) 



a(td 
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tdc 
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(10) 



where tdcc — r is the time when S decays, t dom ~ 
rn p\( av ) 2 i s the time at which it dominates the energy 
density of the universe, and if. , is the time at which it 
freezes-out. Substituting gives 



«(/) = 



1/2 



m-^T-^(av)-^a(t ) 



(11) 



Using this result and Eq. © we again obtain Eq. (JSJ. 

The above discussion is approximate and requires that 
S completely dominates the energy density of the uni- 
verse. Obtaining the perturbations when S does not 
dominate requires that we include the matter contribu- 
tion to the scale factor or energy density during radiation 
domination. This is done in Appendix 1X1 using a differ- 
ent formalism. In Appendix El we confirm these analytic 
results using a numerical calculation of the density per- 
turbations using Boltzmann equations. 



III. EXPLICIT MODELS FOR COUPLING 5" TO 
RADIATION 

It is important to verify that models exist which 
exhibit the features discussed in the previous sec- 
tion. We present two models in which the annihilation 
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cross section is determined by renormalizable and non- 
renormalizable operators, respectively. 

The first model is given by the Lagrangian 



C 



-S 2 - 



(12) 



We assume that X is in thermal equilibrium with the re- 
maining radiation and that S particles are only produced 
through their coupling to X . The interaction terms in the 
above model yield an S decay rate and cross section 



where 



r 



M 



9 m s , 



A 2 
M 2 



T72> (13) 



T when T > ms 
ms when T < ms 



(14) 



Note that we neglect the 0(g 4 ) contribution to the cross 
section. This is justified given the limits on the coupling 
constants derived below. 

The requirement that T t hcrm > ms and that S remains 
in thermal equilibrium down to T ~ ms gives the condi- 
tion on the coupling A 



A > 



TOpl 



On the other hand the condition Tdoc < Tdo 

,3 



9 2 X 4 < 



id'; 
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(15) 
implies 

(16) 



Thus a necessary (but not sufficient) condition on g to 
satisfy both Eq. (JT^J and Eq. (JTSJ is 



.9 < 



m s 

TOpl 



(17) 



Finally, we require that the period of S domination does 
not disrupt big bang nucleosynthesis (BBN). Thus the 
decay of S must reheat the universe to a temperature 
T T h > T B bn, where T rh ~ y/Tm p i. This gives 



T 2 

g2 > 1 BBN 

m s m v \ 



(18) 



Using Tbbn — 10 — TOpi, the above relations provide the 
constraint mg > 10 _21 rn p i. Given any ms satisfying this 
constraint, limits on A and g are calculated using Eq. 1(15(1 
and Eq. JTHJ. 

Note that in this model the S particles are produced 
at T = Tthcrm and remain in thermal equilibrium with 
the radiation until they freeze-out at T ~ m<j. This 
is different from the assumption made in where S 
starts in thermal equilibrium and decouples while still 



relativistic. In order to achieve this scenario, the coupling 
of S to radiation has to proceed via a higher dimensional 
operator, or in other words via the propagation of an 
intermediate particle with mass much greater than ms- 
This brings us to our second model. Consider a heavy 
fermion ips and a light fermion ipxi coupled via an addi- 
tional heavy scalar <\>n with mass mn, 



Ant = gs ^s^s 4>h + gx i>xi>x 4>h 



(19) 



We also assume that the fermion ips decays to radiation 
with rate T. The annihilation cross section is given by 



(20) 



where M is defined in Eq. 114(1 . In this case, thermaliza- 
tion occurs for temperatures bounded by 



m s 



' H 



to|to p i g 2 s g 2 x 



1/3 



< T < g 2 s g x m p i 



(21) 



The conditions that S is in in thermal equilibrium when 
it reaches T ~ ms gives the condition 



g s g x > -3^- 

m s m p \ 



(22) 



Note that one still needs to have a decay rate that is 
small enough such that ^5 decays after it dominates the 
universe. The point of this second example is to show 
that in non-renormalizable models the heavy species can 
either decouple while non-relativistic or while relativistic, 
depending on whether Eq. (|22H is satisfied or not. 



IV. MODELS FOR PRODUCING THE 
FLUCTUATIONS 

The density perturbations in the DGZK mechanism 
and our generalization originate in fluctuations in a light 
scalar field \- I n this section we write down explicit mod- 
els for couplings between S and x- The reason for doing 
this is that these interactions can give rise to back reac- 
tions which can constrain the magnitude of the produced 
density perturbations. Similar results hold for couplings 
between ips and %. 

We find it convenient to define 5\ = X ~ (x)- Note 
that this does not correspond to a perturbative expan- 
sion. The fluctuations in x are created during the in- 
flationary era with 5x ~ Hi n {. Then the leading order 
equation of motion for x can be split into homogeneous 
and inhomogeneous parts, 

(X) - -3H{x) - (V) , 

S X = -3H5 x + H(x)-SV -2<t>(V). (23) 

Here SV = V'— (V), where V is the potential of x and 
the prime denotes a derivative with respect to x- Also, <p 
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is the time perturbation in conformal Newtonian gauge. 
The terms proportional to <f> enter into the leading order 
equation of motion for Sx because their homogeneous 
coefficients do not. 

To simplify the analysis, we first consider the scenario 
where (x) is negligible. From Eqs. I)2i'>ll we see this is 
the case when (%} < 8\- Thus we require the equation 
of motion for (x) to be Hubble friction dominated for 
(x) < This gives the condition 



HH X > H 2 ( X ) > (V). 



(24) 



The fluctuations 5x persist so long as the equation of 
motion for Sx is Hubble friction dominated. With (x) < 
$x this translates into the condition 



H 2 S X > 6V' + 2<f>{V'). 



(25) 



Note that we can combine our simplifying condition that 
(x) be negligible, Eq. (|24|l . with the condition that the 
fluctuations in Sx be Hubble friction dominated, Eq. (|25H . 
Adding these two equations and dropping factors of 2 this 
gives the single condition 



H 2 S X > V. 



(26) 



The constraints of Eqs. i|30|) provide the same upper 
bound to both terms in Eq. (|28l) . Thus the back reactions 
of L x limit the level of density perturbations produced 



via this mechanism to 



c 



< 



< 10" 



(31) 



'pi 



where the last limit on H- 1TL i/m v \ is measured by the 
WMAP collaboration 0. 

The fluctuations resulting from the second interaction 
in C x are linear in Sx and are therefore predominantly 
Gaussian in their distribution. Since the observed level 
of Gaussian fluctuations sets ( ~ 10 -5 , this interaction 
cannot provide a significant fraction of the observed den- 
sity perturbations. However, the fluctuations resulting 
from the first term in C x are quadratic in Sy and there- 
fore non-Gaussian ^(j- Recent analysis 01 limits the 
amplitude of non-Gaussian perturbations to about 10~ 8 . 
Thus we see our model can provide non-Gaussian pertur- 
bations right at the limit of current observation. A lower 
level of perturbations is obtained by reducing as or /is. 

As a variant on the above scenario, we next consider 
the non-renormalizable couplings 



We consider the constraints this condition imposes on 
models for transferring x fluctuations to the radiation. 
We first consider the renormalizable interactions 



■ T S X - T S X 



(27) 



AXL S 2 X 2 -^-?-S 2 X 2 
4 Ml 4 M 2 



(32) 



When x fluctuates these interactions result in fluctua- 
tions in (av) 



and neglect any couplings between x an d A as they are 
irrelevant to our mechanism. When x fluctuates these 
interactions result in S mass fluctuations of 



S rn 



asSx 2 , Ms^X 



4m| 



2m% 



'fa s H 2 



inf 



( US Hint 
V m S 



(28) 
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M 2 



2Sx 
Mi 



K Ml 



Mi 



(33) 



As above, we require that Eq. I|26|l be satisfied. For 
the interactions of C' x this gives 



H 2 S X 



> 



A /2Sx J_ 
4 \M 2 + M 2 



(S 2 X 2 



(34) 



where in the second line we have estimated the size of 
the rms fluctuation at two widely separated co-moving 
points. This mass fluctuation gives rise to fluctuations 
in the decay rate and the annihilation cross section of S 
according to the mass dependence of Eqs. (|13|l . 

As described above, for this fluctuation to persist and 
for (x) to remain negligible requires that H 2 Sx > V . 
Although we assume the self interaction of x is always 
negligible, the interactions of C x contribute to V and 
provide the constraint 



As in the previous example, this constraint is tightest at 



T = m s when {S X )/H 2 



m 



pi- 



Therefore we find 



Ml 



< 



1 1 
Xm 2 pl 



1 

~M 2 



< 



1 H; 



inf 



A m 2 pl 



(35) 



Analogous to the previous example, the constraints of 
Eqs. I|35(l provide the same upper bound to both terms 
in Eq. (|33|l . Thus the back reactions of C' x limit the level 
of density perturbations produced via this mechanism to 



H 2 S X > 



(?*x + f)<**> 



where (S 2 ) is evaluated in the thermal bath, 
straint is tightest at T = ms when (S 2 )/H 2 <■ 
Thus we obtain the constraints 



a s < 



m| 



Ms < 



-Hi 



inf 



'pl 



(29) 



This con- 



m; 



(30) 



c 



>(<jv) 



< 



A m 2 pl 



(36) 



This bound is significantly weaker than the bound of 
Eq. (|31l) obtained via a fluctuating S mass. For exam- 
ple, the fluctuations resulting from C' x could form the 
dominant contribution to the observed density perturba- 
tions if A is sufficiently small. In addition, for a given £ 



5 



decreasing A allows for a lower scale of inflation. Con- 
straints on the smallness of A arc discussed in Section 
PT1 Of course, a lower level of Gaussian (non-Gaussian) 
perturbations is obtained by increasing M2 {Mi). 

Above we have taken (%) to be negligible, which corre- 
sponds to taking (x) < S\- Although this simplifies the 
presentation, it unnecessarily strengthens the constraints 
on /ng and M 2 . We know from Eq. that 



(X>> 



H 2 



(37) 



therefore keeping (x) small implies constraints on the 
potential V. Referring to the second of Eqs. 123(1 . we see 
that for arbitrary (x) the requirement that Sx remains 
Hubble friction dominated gives 

H 2 S X > <P(x) , H 2 5 X > SV , H 2 S X > 4>{V) . (38) 

The first condition provides the constraint (x) < 8x/4>i 
with the evolution of cj> ~ C, described in Appendix ^ 
It is sufficient to take 4> ~ 10 -5 , which also ensures that 
the homogeneous correction that (x) provides to does 
not change raj by more than order unity 2 . Through an 
analysis analogous to that above, we find the conditions 
of Eqs. (|38() constrain the level of Gaussian fluctuations 
for the respective interactions of C x and C' x to 



Cg < 



Sx fn 2 



c s < 



1 (X) H? 



pl 



A Sx TOpi 



(x) 



< 



(39) 



The additional factor of (x) / Sx significantly weakens 
both bounds on Gaussian perturbations. This allows for 
greater freedom in choosing fxs, M 2 , A, and/or H[ n f. 

Non-Gaussian perturbations originate from the cou- 
plings quadratic in \, Taking these into account the 
fluctuation resulting from C x becomes 



C ~ S m = 



asSx 2 , a s (x)Sx . UsSx 



4m| 



2m| 



2m| 



(40) 



Note that the quadratic term ~ as also gives rise to 
a Gaussian contribution to £. Thus, the non-Gaussian 
perturbations obey the relation 



Cng 



as Sx 



Ms + a s (x) 



Cg 



(41) 



Note that taking as — > the non-Gaussian fluctuations 
can be made arbitrarily small. However, even if the as 



2 In Appendix [X] we find that after freeze-out cf> evolves as <j> ~ 
[ps/p)C{> where Q ~ 10~ 5 is the final curvature perturbation. 
Thus if we consider the scenario where \ fluctuations are trans- 
ferred at freeze-out and x subsequently decays, we may take (\) 
to be constrained by <j>~ 1 at freeze-out, which considerable weak- 
ens the bounds in Eqs. 1391 . However, in this case (\) provides a 
homogeneous adjustment to mj which may be much larger than 
mg. This effect could then significantly alter the constraints 
calculated in Section If If I 



term in the potential dominates, the non-Gaussian fluc- 
tuations are always limited by 



Cng < 



g inf 
"*pl 



(42) 



as can be seen by combining Eqs. I|39|l and (|41|) . As men- 
tioned before, WMAP sets the limit H^/m^ < 10" 8 3, 
thus the non-Gaussian fluctuations are at or below the 
current limits from WMAP 01 . In addition, the ob- 
served Gaussian fluctuations can be produced by choos- 
ing Sx/ (x) appropriately. The second model with La- 
grangian given in Eq. (|32|l is less constrained since the 
factor of A weakens the constraint on Sx/(x)- 



V. CONCLUSIONS 

In [{| it was shown that fluctuations in the mass and 
the decay rate of a heavy particle S, which at some point 
dominates the energy density of the universe, lead to adi- 
abatic density perturbations. In this scenario it was as- 
sumed that the heavy particle decouples from radiation 
while it is still relativistic. 

In this work we have shown that if the heavy parti- 
cle remains in thermal equilibrium until it becomes non- 
relativistic, fluctuations in the annihilation cross section 
of this particle with radiation lead to additional sources 
of perturbations. We have presented two simple toy mod- 
els illustrating this effect. These additional fluctuations 
are generic, unless the annihilation cross section is medi- 
ated by an additional particle with mass exceeding ms- 
If the S particle is stable, for example if S is dark matter, 
then the resulting perturbations are non-adiabatic. 

A simple analytical calculation determines the size of 
the density perturbations from fluctuations in the mass, 
decay rate and annihilation cross section. The fluctua- 
tions due to variations in the annihilation cross section 
are shown to be of similar size as the ones generated 
from the original DGZK mechanism. These results are 
checked numerically using Boltzmann equations in con- 
formal Newtonian gauge in Appendix 151 
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APPENDIX A: EVOLUTION OF DENSITY 
PERTURBATIONS 

In this appendix we determine the evolution of density 
perturbations generated by a fluctuating cross section or 
mass during freeze-out. Unlike other analytic derivations 
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given elsewhere in this paper, the one provided here al- 
lows us to easily follow the growth of the non-adiabatic 
perturbation during the radiation dominated era after 
freeze-out. We work in conformal Newtonian gauge with 
negligible anisotropic stress and use the line element 

ds 2 = a 2 [-(1 + 2<j))dn 2 + (1 - 2ip)5 ij dx i dx : >] , (Al) 

with if) = 0. We track the evolution of perturbations 
using the gauge invariant entropy perturbation S and 
curvature perturbation 

We will assume that the scattering and annihilation 
interactions between S and the radiation conserve total 
particle number. This allows us to obtain a first integral 
of the Boltzmann equations. Conservation of total par- 
ticle number in a fixed co-moving volume implies for the 
entropy perturbation 



S = 5 S 



(S R - 30) 1 + 



ns 



nsa 3 



. (A2) 



where 5 s = 5ns /ns and 5r = Snn/nR are the perturba- 
tions in the number densities of S and radiation, respec- 
tively. Here Ao is an integration constant which vanishes 
in the absence of initial adiabatic perturbations. 

Eq. (|A2|) has a few salient features that we now discuss. 
First note that it admits an adiabatic solution 5 = 
whenever both nji/ns and nsa 3 are constant. This so- 
lution is the familiar 5s — 5r = 30 + constant, with the 
constant fixed by Ao- During freeze-out, however, the 
conditions described above are not satisfied, and entropy 
perturbations S L °' are generated at that time. Using 
Eq. tlj we find 



s t.o. 



(A3) 



After freeze-out the heavy particle no longer interacts 
with the radiation and therefore both 5s and 5r obey 
the perturbed Einstein field equation 5s, r = 30. This is 
trivially integrated to give 



30- 



30 



f.o. 



(A4) 



Thus, after freeze-out the entropy perturbation remains 
constant S — <S f, °\ 

To derive the evolution of the curvature perturbation 
C, we use Eq. Q) together with 5ps = ps{5s + 5 m ) and 
5 Pr = {^/^>)pr5r to find after several lines of algebra 



Ps{S io +5 m ) 
3ps + 4p R 



+ I (4° ~ 30 f o ) . (A5) 



Using Eq. JA2J, the second term on the right is sup- 
pressed by {ns/nR){ApR/ ps + 3) compared to the first 
term and is subsequently neglected. We see that in the 
radiation dominated era the curvature perturbations pro- 
duced during freeze-out are suppressed relative to the 
entropy perturbations. 

Once the the heavy particle starts to dominate the 
energy density of the universe, we find 



(A6) 



in agreement with the results obtained in the main body 
of this paper. If S is the dark matter, such a large ratio of 
entropy to curvature perturbations is disfavored by data, 
which requires \S/(\ < 1/3 0. If instead S decays, the 
entropy perturbations in S are transferred to radiation. 
The resulting curvature perturbation is then 



c 



~:5 



(av) 



6 



(A7) 



Here we have finally included the effect of inhomogeneous 
S decay, which are formally calculated in |j| Q ■ 



APPENDIX B: NUMERICAL RESULTS 

We now solve for the evolution of density perturba- 
tions during freeze-out using Boltzmann equations (see 
for example 0, 0])- As above, we work in conformal 
Newtonian gauge and consider only super-horizon per- 
turbations, neglecting all spatial gradients next to con- 
formal time derivatives. Since we assume that S is non- 
relativistic, the distribution function for S is 



fs = e^ T f c , 



(Bl) 



where / cq is the Maxwell-Boltzmann equilibrium distri- 
bution function. Then integrating the background Boltz- 
mann equations over all of phase space gives 



ns 



- 3Hns 
3HnR 



a(av) n„ 



aiav) n„ 



aTns 
+ aTns ■ 



(B2) 



"cq ll S) 

where dots denote derivatives with respect to conformal 
time, Ji = a/ a and n cq is the equilibrium number density 
of S particles. To describe the evolution of the scale 
factor we use Einstein's field equation, 

n 2 = ^(pr + ps)- (B3) 



3m*, 



To derive the subleading order Boltzmann equations 
we allow for fluctuations in mj, (av) and T as defined 
in Eqs. (JSJ). Note that the resulting fluctuations in the 
number density of the radiation Sr correspond to temper- 
ature fluctuations 5t = Sr/3. Therefore the equilibrium 
number density of S particles acquires fluctuations 

3 



(5 m + 5 T ) + exp 



T 



{5 m -5 T ) -1. (B4) 



When the subleading Boltzmann equations are integrated 
over all phase space we obtain 

,2 

' + 5 ((tv) + 2<5c 



<5s — 30 = a{av) 



Sr - 30 



ns 
a(av)ns 

o 

v n 

a(crv' 



6 



5 S ) 

-aT{cj) + 5 r ) 



(av) 



26 c 



+ a(MJ) ^ (0 + 5 {av) + 25 s - 5 R ) 



riR 



+ aT 



riR 



5 r 



(av) 

5s 



Sr) 



(B5) 
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FIG. 1: Evolution of S and £ in units of 5( av ) as a function of 
log(ms/T). The solid, dotted and dashed curves correspond 
to (av)msm p \ = 10 5 , 10 7 and 10 9 , respectively. The large 
values of {av)msm p \ are necessary to ensure that the heavy 
particle is in thermal equilibrium at To = ms/4. 

Note that in deriving the Boltzmann equations we as- 
sumed 5 cq <C 1. This ceases to be valid once T < 
ms(S m — St)- In this case, however, the factor of nl q 
in front of the terms containing S cq is exponentially sup- 
pressed compared to the remaining terms. 



The independent perturbations in Eqs. (|B5(1 are 6s, 
6 r, and <f>. To describe the evolution of <f> we use the first 
order perturbation to Einstein's field equation, 

a 2 

H0 + H 2 (f> = - (6p R + S PS ) . (B6) 
6m pi 

We solve the above system of equations numerically. 
The effects of a fluctuating decay rate T are well-studied 

so we set T = to simplify our results. For con- 
creteness we also assume S interacts with one out of one 
hundred radiative degrees of freedom and we begin in- 
tegration at To = ms/4 with ns(To) = n cq (To). The 
results for the gauge invariant quantities S and £ are 
shown in Fig. ^for several values of (av) ms rn p \. 

The curvature perturbation £ is negligible compared to 
6/ av \ until the heavy particle contributes significantly to 
the energy density. It asymptotes to a value ( ~ —0.3. 
The entropy perturbation S grows during freeze-out and 
soon thereafter reaches a constant value of S ~ —0.9. 
These results and the other features of Fig. are in good 
agreement with the analytical results given in Eqs. (EH) 
and (|A6p . 
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